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ABSTRACT 

A differential method of orbit improvement utilizing observational 
data is presented for Vinti f s spheroidal solution of the dynamical 
problem of unretarded artificial satellite motion about an oblate planet, 
recently modified so as to permit the exact inclusion of the effects of 
the third zonal harmonic term of the planet’s gravitational potential 
field. The first-order Taylor’s series expansion used for the equations 
of condition is fitted to observational values by an iterated least-squares 
process, producing successively smaller corrections to the orbital 
elements. A mean set of elements, conditioned by the observations, re- 
sults for use in orbital predictions for intermediate time points or for 
later epochs. The differential coefficients in the conditional equations, 
applicable to any type of observational data, are derived analytically 
from the equations of the accurate reference orbit. The method of dif- 
ferential correction may be used for orbits of all inclinations, including 
equatorial and polar, and both a first-order and a lengthier, but more 
exact, second-order treatment for the periodic variables are given. 
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INTRODUCTION 


* 


In order to predict the precise position, at a given time, of a satellite revolving about a planet, 
an extensive mathematical theory of the satellite’s motion and the exact values of certain physical 
parameters (e.g., gravitational constants) are required. If the position of an artificial satellite of 
the Earth is to be determined relative to an observer on the Earth’s surface, then, in addition, the 
accurate geodetic position of the observer is necessary. Any mathematical theory of motion is 
based upon certain constants of the motion which initially must be determined empirically. Once 
approximate values for these constants of the motion are available, then they may be utilized in the 
theory to predict future orbital positions of the satellite. A comparison between the positions pre- 
dicted by the theory and those actually observed will indicate discrepancies, expressed numerically 
as the differences between the observed positions and the respective computed positions. These 
are known as observational residuals, and the fact that they are non-zero is due to several causes. 
One fundamental cause which is, unfortunately, ever-present is the inadequacy of the theory: that 
is, the inability to account mathematically for all the physical forces acting upon the satellite. 
Moreover, there are always errors associated with the observations themselves, because of fluc- 
tuations in the atmospheric density, optical imperfections in the telescope, inaccuracies in the 
reduction of the observational recordings, and the like. Further, the physical parameters required 
by the theory often are not known with sufficient precision, and the geodetic position of the observer 
may likewise be measured inaccurately. Finally, there are errors introduced by the approximate 
values for the constants of the motion used in computing the theoretical position. 

Despite all the inadequacies and short-comings enumerated above, the very existence of the 
discrepancies between theory and observation provides a means of improving the approximate values 
of the constants of the motion. One procedure for improving these constants in an analytic theory 
of motion involves expressing the incremental changes in the positional co-ordinates due to changes 
in the constants of the motion as coefficients having the form of partial derivatives in a Taylor’s 
series expansion. In such a case, the partial derivatives must be determined as explicit functions 
of variables which arise in the analytic theory of motion. This method of orbit improvement is 
known as differential correction. 

The purpose of the present paper is to provide such an orbit improvement method for the 
spheroidal theory of artificial satellite motion. The spheroidal method for satellite orbits, de- 
veloped by Vinti (References 2, 3, 4, 5, and 6), supplies a procedure for calculating an accurate 
reference orbit of any drag-free satellite moving in the gravitational field of an axially symmetric 
oblate planet. In the case of artificial satellites of the Earth, the intermediary reference orbit 
reproduces exactly the zeroth and second zonal harmonic coefficients in the series expansion of 
the geopotential function, and also accounts for more than half of the fourth zonal harmonic. Re- 
cently, Vinti (References 7 and 8) has modified the spheroidal potential so as to allow exact in- 
clusion of the effects of the third zonal harmonic as well, heretofore the major neglected effect in 
the spheroidal orbital theory. (The first zonal harmonic is entirely eliminated, of course, by proper 
choice of co-ordinate origin.) Accounting for the third zonal harmonic directly in the intermediary 
orbit in this manner affords a more accurate treatment (Reference 1) than would be possible 
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through perturbation theory. This paper will present a differential method of orbit improvement 
based upon the modified spheroidal theory that includes the effects of the third zonal harmonic 
term. 
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FUNDAMENTAL EQUATIONS FOR DIFFERENTIAL CORRECTION 

Consider a system of constants of the motion, q £ (i =1,2,..., n ), which are utilized in a 
mathematical theory of motion to predict orbital positions of a satellite. In this context, such 
constants of the motion are generally referred to as orbital elements, and the number, n , contained 
in the system is often six. Denote a positional co-ordinate of the satellite at time t by R(t), using 
subscripts "o" and "c" to distinguish between observed values and values computed by the analytical 
theory using the orbital elements. We assume that the differences of the co-ordinates, R 0 (t) - R c (t), 
as well as the required corrections to the elements in the improvement of the orbit, are sufficiently 
small so that their squares and higher powers may be neglected. We can then express the obser- 
vational residuals by a truncated Taylor f s series expansion restricted to first powers, in the 
general form: 


(*) 

R o(t) - R c (t) - Aq > * 

i=l 

Here the computed value, R c (t) , is viewed as a function of n independent variables, q £ (i = 1, 2, . . . , 
n), which are to be improved by the additive increments, Aq. (i =1,2,..., n). The coefficients 
in the Taylor’s expansion expressing the increment of the co-ordinates caused by a change in the 
orbital elements have the form of partial derivatives, and these must be determined analytically 
from the equations in the mathematical theory of motion. Each separately observed co-ordinate 
yields an equation for corrections of the elements of the form given above. Ordinarily, the number 
of such so-called equations of condition available far exceeds the number, n, of unknowns, Aq. . 

This set of linear simultaneous equations forms an inconsistent system (due to inherent random 
and possibly systematic errors in the observations) for which no exact solution exists. It may be 
solved by the method of least squares, yielding "preferred" values for the unknowns. The solution 
of the equations by the principle of least squares follows well-known schemes (Reference 9), so 
that the problem reduces to evaluating the derivatives of the co-ordinates with respect to the el- 
ements (the so-called differential coefficients). 

We now introduce certain conventions in order to allow discussion of the problem in more 
explicit terms. Assume that the satellite’s positional data are recorded at the tracking stations in 
the form of direction cosines observed with respect to a topocentric (i.e., situated at the Earth’s 
surface) co-ordinate system. The topocentric or "local" co-ordinates shall be distinguished by 
the subscript, "M". The system is orthogonal and right-handed, with the X M - Y M plane tangent to 
the Earth’s surface. The X M -axis extends in an easterly direction along the line of latitude, the 
Y M -axis extends in a northerly direction along the line of longitude, and the Z M -axis is normal to 
the surface and points toward the geodetic zenith. If L 0 and M 0 denote the observed direction cosines 
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of a satellite (in the X M - and Y M -directions, respectively) for a given time of observation, then the 
corresponding computed values of the direction cosines are given in terms of the local co- 
ordinates by 


c = ( V + y m 2 + Z„ 2 ) 1/2 ' 

and 


' = (V + Y M J +Z M 2 ) 1/2 

The computed value of the third direction cosine, N c (in the Z M -direction), is not an independent 
parameter, but is pre-determined by L c and M c through the relation, 


N 


c 


(l-V-Mc 2 ) 


1/2 


For this reason, each satellite observation provides two, and only two, independent co-ordinate 
values, which here are chosen to be L 0 and M 0 . 

The mathematical theory of motion will ordinarily predict a satellite’s position with respect 
to a rectangular geocentric system (x, Y, z) of co-ordinates, the so-called inertial frame of ref- 
erence. In this system, the X - Y plane is the Earth’s equatorial plane, and the origin is situated 
at the Earth’s center of mass. The X-axis extends toward the vernal equinox (the first point of 
Aries), the Y-axis is orthogonally to the east to form a right-handed system, and the z-axis coin- 
cides with the Earth’s polar axis. In order to obtain a satellite’s local co-ordinates from its 
inertial co-ordinates, the inertial co-ordinates of the observation point at the time of observation, 
which we shall denote (x T , Y T , Z T ) , must be known and two rotations performed to bring the local 
and inertial systems into parallel alignment. Let us denote by 4> x the angle between the vernal 
equinox and the X M -axis as measured in the observation latitude plane. Then ip x will depend upon 
the longitude of the tracking station and the hour angle of the vernal equinox at the time of observa- 
tion. If 0 D represents the latitude of the tracking station, then the local co-ordinates of the satel- 
lite at observation time are given by the matrix equation, 


sin# D cos @ D 


cos @ D sin 0 D 


cos 0 x sin \p x 0 


- sin^/ x cos i// x 0 


0 1 


^ 
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The difference of column matrices on the extreme right represents a translation from the Earth's 
center to the tracking station position; the center matrix on the right represents a rotation in the 
latitude plane about the polar axis through an angle, </<„ , to bring the inertial X-axis into coincidence 
with the station's X M -axis; the remaining matrix on the right represents a rotation in the longitude 
plane about the X M -axis through an angle equal to the complement of 0 D to bring the inertial Z-axis 
into coincidence with the station's z M -axis. This matrix equation, when expanded, reads 

X M = (X - X T ) cos \p x + (Y-Y T )sin^ x , 

Y m = - ( X - X T ) sin \p x sin 0 D + ( Y - Y T ) cos \p x sin & D + (z - Z T ) cos & D , 


and 


Z M - ( X - X T ) sin \p x cos 0 D - (Y - Y T ) cos \p x cos 0 D + (Z - Z T ) sin 0^ . 

We can now write the first-order Taylor’s series expansion for the equations of condition 
corresponding to each time of observation in the following more explicit form: 


V" aL c 

l c ^ Aqi 
i=i 


and 


^ M c 

- M c = 2 - • 

i=i 


Since the local co-ordinates are functions of the inertial co-ordinates, which are in turn 
functions of the orbital elements, then the chain rule may be used to expand the differential 
coefficients. 


dh c dh c <?X M dL c dY u SL C (9Z m 

aqT = a^ a^7 + av; a^7 + az„ a^7 

and 

<?m c _ aM c ax M aM c aY M aM c az M 
aqT “ ^^ + aY M a^7 + ^z M aq7 
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The partial derivatives of the direction cosines with respect to the local co-ordinates are found 
directly from the expressions for L c andM c . 


and 


where 


*Lc 

* X 2 

1 a m 

<? X m 

R M R m 3 


X « Y « 


V 

dK 

X «Z« 


r m 3 

<9M C 

X m y « 


r m 3 

aM c 

; J, _ V 

aT u 

r m R m 3 

<9M C 

y m Z m 

IT* 



R m - (v + y « +z « 2 ) ,/2 • 

Since the co-ordinates, X T , Y T , and Z T , and the angles, «/» x and d D , are independent of the orbital 
elements (and merely geodesic functions), we have the matrix equation. 


ax M 


3q. 


5Y m 



II 

dZ„ 


i q L 



0 


sin # D 


- cos & D 


cos 8 t 


cos \p x 


- sin v// 


sin 4> x 


0 


dX 

dq, 


cos \p x 


0 


_dY 

dq, 


0 


1 


dZ 
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Equivalently, 


* 


5X m _ dX dY 

dq _ - cos^ x ^ ' 

3Y m _ . , . „ dX , . a d\ n 8Z 

dq. ' ~ s »n^ sxne p ^7 + cos i/- x sin£> D ^7 + cos 0 D ^7 ' 

and 

aZ M _ . . Q dX dY . a dZ 

dq, " sin^/ x cos e D 1^7 - c °s^ x cos 0 D 5^7 + sine D aqT ' 

The problem remains to evaluate the differential coefficients, dX/dq, , dY/<?q. , and dZ/dq i (i = 1, 2, 

• • • , 6), which are the derivatives of the inertial co-ordinates with respect to the orbital elements. 
At this point, the mathematical theory of motion becomes of primary importance. 


OUTLINE OF METHOD USED TO EVALUATE DIFFERENTIAL COEFFICIENTS 

The constants of the motion, q. (i =1,2,..., 6), which we choose for the mean orbital 
elements for the modified spheroidal theory of satellite motion (including the exact effects of the 
third zonal harmonic coefficient of the oblate planet’s gravitational field) are the following: 

q x = a, the semi-major axis. 

q 2 = e, the eccentricity. 

q 3 = S, corresponding to sin 2 I in Keplerian motion, where I is the inclination of the orbital 
plane to the equator. (However, s may be negative for orbits sufficiently close to equatorial. For 
a polar orbit, S is unity, and in all cases, S 2 < I.) 

q 4 = , corresponding to the negative of the time of passage through pericenter in Keplerian 

motion. 

q 5 a f$ 2 , corresponding to the argument of pericenter in Keplerian motion. 

q 6 = , corresponding to the right ascension of the ascending node in Keplerian motion. 

These elements differ slightly from those selected by Vinti in the final algorithm for the reference 
orbit (Reference 8, Section 12), but the above parameters seem more suitable for the differential 
correction. 
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The relation between rectangular inertial co-ordinates and oblate spheroidal co-ordinates, 
p 9 7) 9 <p 9 which are involved in the solution of the problem of satellite motion, is given by 


X - i{p 2 + c 2 ) ( 1 ~rj 2 ) cos » 
Y - )/( p 2 + c 2 ) (l - t ] 2 ) sin</> , 

and 


Z - prj - 5 , (-1 £17 1 1 ) . 

Here c and 8 are adjustable parameters that are chosen to agree with the coefficients of the zonal 
harmonics in the series expansion of the Earth’s potential function. In terms of the Earth’s equa- 
torial radius, r e , and the n 1 * zonal harmonic coefficient, J n , the proper choices are 

- 2 = ?Js 2 Jr’) . 


and 


S = - J f e Jz" 1 J 3 ' 

For the Earth, the values are approximately c = 210 km and 8=7 km. 

In order to evaluate the derivatives of the inertial co-ordinates with respect to the orbital 
elements, we must know dp/d q i? dv/dv i9 and d<f>/dq i (i = 1, 2, • • • , 6). The oblate spheroidal co- 
ordinates are rather involved functions of the orbital elements. The process for determining the 
partial derivatives of p 9 77, and <f> is a lengthy one which will be presented in a synthetic, rather 
than analytic, manner. That is, necessary partial derivatives of the simpler functions of the orbital 
elements will be given first, followed by partial derivatives of more complicated functions of the 
elements involving the pre-determined partial derivatives. 

It is worth noting in passing that the equations to be presented herein apply to orbits of all 
inclinations, including equatorial. There are no special simplifications introduced in the case of 
equatorial or near-equatorial inclinations, as there were for the spheroidal theory that did not 
include the effects of the third zonal harmonic (Reference 6). 

However, as in the earlier spheroidal theory, the differential correction may include derivatives 
of the periodic terms taken through the second order, or, alternatively, it may be simplified to 
omit periodic terms higher than first order. In either case, some second-order effects are in- 
cluded with the first-order terms and even with the zeroth-order terms. The terms of second and 
higher orders that are to be omitted in the simplified version will be indicated as such. 
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TIME -INDEPENDENT PARTIAL DERIVATIVES IN THE DIFFERENTIAL CORRECTION 


We begin with equation (23.3),* which reads: 


p - a(l -e 2 ) . 


Thus: 


and 


dp 

da 


1 - e 2 , 


dp 

de 


2 a e , 


<?P _ 

as " 


0 . 


Since 


b i 



by line 8 of section 12 of Reference 8, then by equation (40.1), we find: 


ac 2 (ap - c 2 S)(l-S)-4a 2 c 2 S 2 p l S(l-S)[l + c 2 (ap) 1 (3S - 2)] 
k' (ap~c 2 )(ap-c 2 S)+4a 2 c 2 S + 4c 2 S 2 p~ 2 S(l-S)(3ap-4a 2 -c 2 ) 

If we denote the numerator and denominator of b, by N and D respectively, then we have: 

5b, f 

~da = - c 2 (1-S)D-‘|c 2 S 

+ 4 (^) 2 S 2 S(1 -S) |"3p-8a+ (s ~ + 2 ^ - 3 a 


.(2p + a|f) + 48 2 s[j (2-||f) + g) 2 (3S - 2) (l - 2 


a ap ' 
p da , 


ND 


2 2 ap - c 2 - c 2 S)^p + a + 8 ac 2 S 


*A11 equation numbers, used in specifying the defining relation for a given variable, will refer, until otherwise indicated, to 
erence 8. 


dp 

da 


Ref- 


8 


db 1 

de 


: 2 D- 1 (l-S)|f 


1 + 4 (p) 2 s[l + : 


ap 


(3S - 2) 


NIT 2 


de 


>(2 ap - < 


s) 



S 2 S(1-S) (s^+2j-3a) 


and 


abj 

~dS 


- ac 2 D 


2 n-i 


ap + (1 - 2S) (4 J S 2 + c 2 ) - 4(|j S 2 (9S 2 - 10S + 2 ) 


- ND " 2 


[c 2 (4 a 2 - ap + c 2 ) + 4^) S 2 ( 1 - 2S) (3 ap - 4 a 2 



Since 


b 2 - /B , 

by line 8 of section 12 of Reference 8, then by equation (40.2), we find: 

b 2 2 = c 2 - bja" 1 (ap-c 2 ) . 


Thus: 


ab 2 

<9a 


j (*.)•■ 




dp \ 
da) ’ 


and 


db 2 

de 


j. 

2 


V 1 



+ b 


dp 
1 de 


db 2 

~dS 


- \ (abj)' 1 (ap-c 2 ) 


ab t 

"as" 


This last sequence of three equations can be written as a single generalized equation if we intro- 
duce the following notation. Let q 2 = a, q 2 = e, and q 3 = S, and let 8 i . be the Kronecker delta, 
defined as follows: S.. = 1 when i = j, and 8.. = 0 whenever i / j . Then, for i = 1, 2, 3, 


db 2 

~d^[ 


3 ( ab 2 ) _1 



b , a' 


abjN / 

aqTj ~ b l \P S li + a 
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This notation will be used extensively in order to write later sequences of equations in an efficient 
manner. 

Combining equations (35.1) and (35.2) with equations already mentioned, we find: 


Po 


(a + bj ) 1 (b 2 2 + ap + 4 abj - c 2 ) . 


Thus: 


and 


dp 0 

~5Z = ( atb i) _1 


■’Po , > . 

i9e ( a + b i) 


[ db 2 / 

<5b l <?p\ ( 


2 b 2 ~da + P + 4b l +a V 

1 da + <9a j Po \ 

1 da j 


db 2 

2 b 2 ~Te + a 


( ^ b i a P \ 

\4~te + de)-Po Sv J 


^P 0 

Is z ( a+b i)' 1 


db~ 


<9b x 


L 2b 2 dS + ( 4a “Po) <9S 

This may be written in the generalized form, for i = 1, 2, 3, 


3p 0 

n£ = (- +b i) _1 


dh 2 , . / db l <3p\ / 

2 b 2 + ( 4b i + p) s u + a ( 4 + " Po^n 


c?b j 


From equation (21.2), we have: 


Thus, for i =1, 2, 3, 

i 


d a 2 (j, <?P 0 

dq i ~ 2 cl 2 <9q. 


For convenience in later equations, we now differentiate the inverse of 


a oPo = ( a + b i)Po 
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to find, for i = 1, 2, 3, 


<9 

dq\ 


( a oPo)' ! 


- p 0 -‘ (a + bj)-‘ 



From equations (31.1) and (31.2), neglecting terms of fourth order, we find: 


_1 

U 


(l-S) 


(28 . 

r.) “' s > 


i* 


(‘ - Hk s ) 


(1-2S) 


Now, for i = 1, 2 only, we have: 


<ju 

d q, 


- u 2 (l-S) 



d 

<?q, 



c 

i oPo 






a oPo 


(1- 2S) 


2c 2 



(1- 2S) 




^- 2s )] A 


( a oPo) 


However, 


du 

dS 



-2c 2 



(l-S) 




a oPn 


l" 3 

( 1 “ 2S)J [d-2S) ^ (a 0 p 0 )-* -2(a 0 p 0 ) 


By equation (32.1), 


C 


2 


a oPo 


u 


Thus, for i = 1, 2, 3, 


_ d . v . c 2 <9u 

dq t ~ c u dq i ( a oPo) + a 0 p 0 ^q. 
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By combining equations (32.1) and (32.2), we have: 


2b 


(l-C 2 )(l-C 2 S)-‘ 


Thus, for i = 1, 2, 3, 


^7 -- I (i-c 2 s)- 1 {(i-c 2 )(^ -^^t)-u(i-c 2 s)-» 




(1-S) -s 3i (i-c 2 )c 2 


By combining equations (32.1) and (32.3), we have: 


P - — u(l-S) (l-C 2 S)-> ■ 


Thus, for i = 1, 2, 3, 


dP 

dq, 


<9u 


1 ^2^) 1 l ^ 3i ^2 + ^ <9q i J Po dq { +u dqi 


S 3i ( 1 C 2 S ) 1 P 0 U 


Equation (21.1) may be rewritten as 


M(a +b 0 _1 


Thus, for i = 1, 2, 3, 


da 

Jq i 


1 1 / \ 

2^{ a+h i) 


db. 

- 2 i 8 u + ^ 


For convenience in later equations, we now evaluate, for i = 1, 2, 3, 


_i_ (h s 

5q, U 


5b j 5b 2 

b 2 1 dq~ b l b 2 


and 


d ( h 2 

a Qi 


db 2 


, - 2 dp 

" b 2P lq[ 
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By equation (5.31),* we have: 


= (l 


e 2 ) 1/2 p 


LW 

n - 2 


P n (^) R n - 2 [(1 - e 2 )‘ /2 ] 


Here, and in the following, P n (x) is the Legendre polynomial with argument x of degree n. The 
definition of R is given by equation (5.28), viz., R n (x ) - x" P n ( i/x) , where 0 < x < 1. We shall de- 
note by P n ' (x) the derivative of the Legendre polynomial with respect to the argument. Then for 
i = 1, 2, 3, 



1 5 2i P° 




n-3 



n-2 


[ (l _ e 2)-./ 2 ] 



By equation (5.36), 


A 


2 


(l “ e 2 ) 1/2 p 1 


CD 


n = 0 


h\ 

V?J P n ^ ) R n 


[(l -e 2 ) 172 ] ’ 


•Equation numbers used for the defining relations will now refer, until otherwise indicated, to Reference 4. 
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so that, for i = 1, 2, 3, 


d = - A 2 P'‘ ' s 2i Ajetl-e^-'+fl-e 2 ) 172 ?- 1 


5 ( b 2\ f- /b 2 \ n_1 /b, 

5q i \ P ) Z_ l n \ P / P " l b 2 


n El K [(l-e 2 ) 1/2 ] 


d ( h l\ t^( h 2\\,( h 


dq, (b 


+ 8,, p 1 e 


n'(bl R n [(l-e 2 ) 1 


{(1-e 2 )- 1 £](t) [(l- e2 ) 1/2 ] n P n (^)p„' [(l-e 2 )' 

L n= 1 


£"(#) [d-e 2 ) 1 / 2 ] nM P n (^)p n [(l-e 2 )-‘ /2 ] 


By equations (5.61), (5.50), and (5.53), we have: 


A 3 = (l-e’)*"p-» £ D n R n+2 [(l - e 2 ) 1 


where D is computed as follows: 


D „ = °2j 


r-» /c\ 2i 2m / ^ 2 \ m 

/ , ^ _1)) m (p) \P/ P 2m 

m = 0 


(n an even integer), and 


D n = D 2j +1 


£ (I) 


2j-2m /b 2 \ 2m + 1 / bjN 


p J r 2m + l l b 


(n an odd integer). Then, for i = 1, 2, 3, 


A 3 3P' 1 W. +S 2i e 


03 

(l - e 2 ) " 1 J ♦ (l-e^^p-3 £* n+2 [(l-e 2 ) ,/2 ]^7 


ep" 3 A ( 1 - e 2 ) -1 2^[(l-e 2 ) 1/2 ] n+2 D n P„; 2 [(l-e») ! 


£](n + 2) [(l ~ e 2 ) 1/2 J " +1 D n P n+2 [(l-e*)' 
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where <9D n /^q. is computed as follows. If n is an even integer, then 


dq, <9q i 


2cp 


^ (f) 


2(j-m)-l / b ^ 2m 


P 2m V b 2 , 


+ 2 ^_ (M f,.nr. f b i 


5qT 


2m\ b 


5 / b i\ v -1 ^ b N2m 

W L ( " 1), ' ra fe) 


p ' 

P / 2m 


If n is an odd integer, then 


<?D„ <9D 


2 j+i 


^ Q i <9q. 


n - 2 

' 2C P ^ 


i£<- 


-m) o 


2(j-m)-l /b 2 \ 2m+1 


P 2m + 1 \bjJ 


JL ( h _A - 

*>i \P 


\ 2 (j “m ) ^ b 3 \ 2m /b, 

\ P / I- 2m+l\b, 


R) P. 


_a_ /_V J 

\b 2 




2(j-a.)/b 2 \ 2m+1 


P ’ 

2m + l 


By equation (5.37), we have: 


A 21 l 1 e 2 ) 172 p 2 e|b, + (3b, 2 b 2 2 ) p 1 2 b l b 2 2 p 2 ( J + 4 e2 ) + b 2 4 P" 3 ( 4 + 3e 2 ) 

Thus, for i = 1, 2, 3, 

~aq“ = _ A 21 Ip " 1 ^ + S 2i [ e (* -e 2 )' 1 -e' 1 ]^ (l -e 2 ) 1 / 2 p“ 

(9o / ^b, db 2 \ 

- 2p- 2 (3b, 2 -b 2 2 ) +2p-> (3b, -b 2 ^ 


k -1 , dbl 

~ b l P ^ + 


+ |b 2 p- 2 


- dp ^ D i do 2 

b 2 P dq, ~ b 2 <9q7 " 2b l dq7 


(l + l e 2 )( 3 b, 

b/p -3 ( 4+3 e 2 )^ -b.p’ 1 ^7) + § 8 2i b 2 2 p- 2 e(b ; 


3 

+ 7 


, 2 P M -b 
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By equation (5.38), we have: 


A 22 - ■§" (l " e 2 ) 172 p ~ 3 e 2 [(3bj 2 - b/) ~ 9 bj b 2 2 p 1 + ^ b 2 4 p 2 (b + e 2 )] 

Thus, for i = 1, 2, 3, 


dA 


22 

^q“ 


-i dp 

P +S 2. 


I . 

4 b 2 i 


> 3 e 2 

( 1 - e 2 

) 1/2 j 

[[-»■ 

1 (3b,' 

dp 

+ (ab, 

9 

_1 , » 

\ 5b : 

dqT 

- 2 P 

' V 

) 


>b,b 2 

_ 

1 


- 9p" 

” b 2 _ 

1 <59 7 ] 


e (l - 

e 2 ) 1/2 

[3b , 2 

-b 2 2 - 

9p ~ 1 
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By equation (5.62), we have: 

*3, -- (l-e 2 ) 1/2 p- 

Thus, for i = 1, 2, 3, 


>i p ' 1 (3 


2 + b, p ” 1 13+^-e 2 ) - p 2 (-5 


2 b 2 2 +c 2 


) (4 + 3e 2 ) 


<?A 


31 

dq. 


- - 3 A,, p 


1 


3 lP <9q ; 


(l - e 2 ) 1/2 p 4 e " 


{ “bj P ' 1 (3 + | e 2 ) 


2P ' 2 ( 4 + 3e 2 )(-jb 2 2 + c 2 ) + (. 


. 3 ,\ dh i 

l 3+ 4 e i ^ 


- S„ 


db 2 i /1 

- b 2 p ' 1 (4 + 3 e 2 ) 3^7 + S 2 . |j b, e - 6 p ~ 1 e(^ b 2 2 + c 2 


j p 3 Jjl " e 2 ) 1/2 e 2 - (l - e 2 ) 1/2 ^J 2 + b, p 1 (3 + e 2 ) - p -2 (-5 b 2 2 + c 2 ) (4 + 3e 2 ) 


By the corrected version of equation (5.63), 

[? + i b j p ' 1 -p ' 2 (•§ b 2 2 + c 2 )(| + ^ e 2 )_ 


A 32 = ( 1 - e 2 ) 


1/2 p 3 e 2 
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Thus, for i = 1, 2, 3, 


«9A 


3 2 


( 9 q. 


- A 


32 


3P_1 + S 2i e ( X _e2 )" 


(l - e 2 ) 1/2 p 4 e 2 


P' 2 (l b 2 2 + c2 ) ( 3 + 3>e 2 ) 


3 , 

4 P b i 


_3 P 3 _1 _ 

+ 4 3q ; " 2 P b : 


( 3+ 5 e2 ) 


3 b 2 


S 2i (l- e q^p- 3 e [| (l + 3 bj p“ 1 ) - p ' 2 


By the corrected version of equation (5.64), 


‘•33 = 3 (l- e2 ) 1/2 p- 4 e 3 


l b l “P" 1 [\ b 2 2 +r2 ) 


Thus, for i = 1, 2, 3, 


3A„ i 

~3q7 = ” 3A 33 P" 1 3q: + 3 ( X _e2 ) V2 P' 4 e3 


2 P 2 (l b 2 2 + c2 ) " ? P"‘ b i 




1 ab l 

+ 4 - P b. 


dtO i r 

3q7j + 3 S 2 i P' 4 e ' [ 3 ( X _e2 ) 


1/2 


- (l-e*)-va e 2 ] 


1 

4 b i 


By equation (5.65), 


A 34 ' ~ 3T ( 1_e2 ) 1/2e4 P" 5 (l b 2 2 + c2 ) 


Thus, for i = 1, 2, 3, 


thl _ -! ip 1 /, 2 x 

dq { ~ “ 5A 34P 3qT ~ 32 l 1_e ) 


3b, 

1/2 P ~ 5e4b 2 3q~ 


f S 2i e3 P’ 5 [k (l-e’)-" 2 e’-(l-, 


b 2 2 + c 2 ) (3 + e 2 )^ • 


-p ’ 1 (j b 2 2 + c 2 ) 


2 H(K 2 + c 2 ) • 
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By equation (47),* we have: 


Q = y' P 2 + S . 


Thus, for i = 1, 2, 3, 


_5Q 



1 

+ 2 


8 


3i 


)• 


From equation (102.1), 


l = P(l-S)-> 


so that, for i = 1, 2, 3, 


ai 


(i -s )- 



s 3i p (i- s ) 


From equations (147.1) and (147.2), 

h , = J ( i + c ^)-" 2 


and 


5 U" c i -C 2 )- ,/2 • 


Thus, for i = 1, 2, 3, 


3/^1 , 

5qj 2h l \dq ; dqj 

and 

ah 2 /ac, ac 2 \ 

dq~ ~ 2h 2 3 \dq7 + 5q7/ 


* Equation numbers will now refer, until otherwise indicated, to Reference 8. 
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From equations (100.1) and (100.2), 


Thus, for i = 1, 2, 3, 


e 2 = O(l-P)' 1 


e 3 = Q(1+P)' 


-l(J± J>p 


= (1-P > +e 2 


<9e 3 _ _./<JQ dP\ 

” ( X tP > (<?q. e 3 <9q i j 


From equation (131.2), 


so that, for i = 1, 2, 3, 


: ( a+b i ) _1 




+ s 2i a-e' IS j ; +7q: 


From equation (149.1), 


a 3 = (sgn a 3 )a 2 (l - Su' 1 ) 1 


where sgna 3 = a 3 j ja 3 | is +1 for a direct orbit and -1 for a retrograde orbit. Then, for i = 1, 2, 3, 


d a 0 1 / 

sgn a 3 (l - Su _ 1 ) 1/2 ^ a 2 u ' 1 (l - Su“ 1 )" I/2 ^Su“ 


< — -8 
dq. b 3i 
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By equation (76), neglecting terms of fourth order, we have: 


B/ = \ Q 2 + P 2 - J C, PQ 2 + | C 2 P 2 Q 2 + («2 + 3C i 2 ) Q 4 


t! C i C 2 PQ 4 + 256 (2C 2 2 + SCi 2 C 2 ) Q 6 + 2MS C 2 3 Q 


15 


175 


Thus, for i = 1, 2, 3, 


^ - q[i + Q 2 (4C 2 + 3C 1 2 ) + tS c 2 Q 4 (2c 2 + sc 1 2 )]^ + 2P^ -Iq^pq^: +c 1 Q^ 1 i 


<?Q 


dP 3 


a *i 


aCj 


dP 


_dQ\ 3 

"1 r c? <q i 


dC n 


dP 


dQ 


dC , 


ac, 


+ 2 C, p I + |pq!pQ ^7 + 2 C 2 Q^- + 2 c 2 p^j + 31 Q 4 ^ai; +3C i a q , 


45 ( dC ac 2 a p dQ 

3 § Q 3 (c 2 PQ 5^7 +C lP Q a — +C 1 C 2 Q a - Ti + 4C, C 2 P ^ 


15 

+ 256 Q 


ac 2 _ ac, 


|^(4C 2 + 5Ci 2 ) ^q. + 10Cj C 2 ^q. 
By equation (65), neglecting terms of fifth order, we have: 

B 2 = 1 -1 C,P + (| Cl 2 + ^C 2 )(p 2 + ^Q 2 ) + ^C/Q^ -| Cl C 2 PQ 2 


175 

2048 


/ <?C 2 dQ \ 

C/Q 7 ( 3 Q d~q[ + 8C 2 ^J 


"" 32 c i 3 PQ 2 + TT8 C 1 2 c 2 Q 4 + ^ c 2 3 Q 6 + S C 1 4 Q 4 


+ ¥ c 2 2 P 2 Q 2 - T2§ C 1 c 2 2 PQ 4 + C 1 2 c 2 2 <? 6 + (m ) 2 c 2 4 Q 8 

Thus, for i = 1, 2, 3, 


+ 32 C 2 Q3 
ac 


75 , 1 

/ 35\ 

2 i 

/ ^ 2 


9 Jr 

[9 + g-C 2 Q 2 + 2 

\ 8 ) 

c 2 2 Q 4 

V Q 

+ ‘ 2C 2 3q, J 

" 8 Q \ ( 


ac. 


+ C .PQ^ + C lC2 Q^ + 2C,C 2 p|) - H Cj 2 Q ^3PQ ^ + C, Q + 2C, P ^ 


45 


ac t 


ac. 


aQ \ 105 


128 C, Q 3 (X 2 Q ^ + C, Q 357 + 4C, C 2 + -jgg C i 3 Q 3 (Q 5q, + C i 3^ j + ¥ C 2 ^ a^T 


ac, 


aQ 


ac. 


ap 


aQ\ 225 


c 2 Q aqT +C 2 p acfj - T28 c 2 Q 3 ^c 2 pq + 2C,PQ^ 


ac. 


ac. 


c,c 2 q^ + 4 Ci c 2 p^) + +c lQ ]^ + 3c,c 2 g 


dC , 


dC , 


0 
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By the corrected version of equation (95), neglecting terms of fifth order, we have: 


1 

3 , 

/ 35 , 

15 , 3 ,\ 

( 1 9) 

3 

B 3 = ” 2 C 2 

- 8 C i 2 - 

Via c i 4 + 

16 C 1 C 2 + 8 C 2 ) 

U + 2 Q 2 } 

8 


- Ci c 2 + K c »)(' + W+iv)- ^ c 2 4 ( i+ Iq 2 + Iq 4 + ^q 6 ) 

+ {re C ! 3+ 1 c i c 2 ) P + if ^/(p+IpQ 2 

Thus, for i = 1, 2, 3, 

$ ■ C | ^)-( 1 ^Q>)[! c ,Sc ? *3C 2 )^.5(| CiltCi )^] 


2 Wq ; + 2 C l 


<v + il^c, +4^-41^^^ + 4 Ic/(i + 4q 2 + ¥q 4 )Jq^-I c 2 p ( p ^: +c 2^ 


TS c 2 li + 3 Q 2 + I Q 4 ) ft c i c 2 + (“2 c i 2 +c 2 ) Sri + 4 c 2 (1 + 1 Q 2 ) |3p(c 2 ^ 



dP 

(21 , „ > 

dQ 

35 - / 

2C ' ^r.l 

+ 3 C C ^ — 
^2 dq i 

- c 2 q[ t c , 2 + c 2j 

^ ^qi 

32 C 2 l 


1 , 3 . 5 ,\ ^2 

1 + 2 Q + 8 ^ + 16 Q ) <9q ; 




By equation (116.3), 


i Q 2 - I C 2 Q 4 


so that, for i = 1, 2, 3, 


= -^Q(i +p 2Q 2 )S-^ 4 S • 


-5^7 - - 2 '^2 v 

By the corrected version of equations (116.4), 


B 21 = - C 2 P 0 + ^ CjC.Q 3 


21 


and 


B 22 - -12 [l«2 + 3C, 2 )Q 2 + 3C 2 2 Q 4 J . 


Thus, for i = 1, 2, 3, 

i9B 21 1 

dq, - 2 

and 


9 -> 


/ 9 „ \ 

X 2 

/ 27 

, 1 ) 

+ 8 C 2 Q 2 y 

^ 5q; 

+ Q\i6 C 1 Q 2 - p ) 

5qT + 

lie C i 

c 2 Q c 2 p + j Cj ) 


IQ . c o i* 
aqi L 2 ^ ( 9 q i 


<9B 


22 


d *i 




dC, 


SC, 


f90 t/v '2 uv i 

(4C 2 + X, 2 + 6C 2 2 Q 2 ) + Q (2 ♦ 3 C 2 Q 2 ) ^ + X, Q ^7 

By equations (122.1) and (122.2), we have: 

2nv l - (-2 ai ) 1/2 (a +b, + A, +c 2 , 

and 

277V 2 = a 2 u-‘/ 2 A 2 B 2 -‘ (a + b, tA, *c 2 A 2 B 1 'B 2 - , )-> 

Thus, for i = 1, 2, 3, 

277v j j(a + bj + A. + c 2 A 2 Bj' B 2 1 ) 


W- (^ 1 ) 


dbj 

’*> + 


dk 2 t9B,' SB 

+ Sq~ + ° 2 B 2 * ^ B i <5q7 + A 2 c)q~ ” A 2 B l B 2 dq ; 


da. 

- n 1 

2 a i dq i 


and 


~d^. (^ 2 ) = - 2t7v 2 j(a +bj +Aj +c 2 A 2 B/ B 


db 1 d\ f <9A 2 

5 ii + + + c 2 b 2 - 1 (b/ ^7 


<9B/ dB 2 

+ A 2 ~^q x ~ A 2 ®i B 2 1 <9q i 


doL„ 


dA 


2 ‘ dq A + A 2 ‘ <9q £ 2 u " dq A c 2 


2 1 _ , du 


dB, 


- EL 


22 



The following time-independent partial derivatives are used only when the differential cor- 
rection includes periodic terms through the second order. 

By equations (5.32), (5.33), (5.39), and (5.40) of Reference 4, we have: 


A 11 " 4 ( X “ e2 ) 1/2 P" 3 b 2 2 e ( b 2 2 " 2b l p) 


1 2 


- 32 (l - e 2 ) 1/2 p" 3 b 2 4 e 2 , 


A 23 = (l - e 2 ) 1/2 p” 4 b 2 2 e 3 (b 2 2 p -1 - b,) 


and 


A 24 = 2§§ (l _ e 2 ) 1/2 p" s bj 4 e 4 


Thus, for i = 1, 2, 3, 


r)A r 

= “ A u { 3 P _1 Tq ; + S 2. [e(l-e 2 )-' -e’ 1 ]) 


T d n ^ b t 5b a "| 

(l _ e 2 ) 1/2 p‘ 3 b 2 e[b 1 b 2 + 2(b, p - b 2 2 ) 


<?A 

~3q 

dA, 

~dq 


if = I (1 " e2 )‘ /2 P' 3 »>/ e2 (Sf - I P“ b 2 + A *2i P* 3 eb/ [(1 ~e 2 )" 2 - i (1 - e 2 )-^ e 2 ] 


7 = " A 23 P' 1 + 1 (l-e 2 ) 1/2 p- 4 b 2 e3 j^p' 1 (3b, - 4 b 2 2 p _1 ) ^ 

(9bi <9b 2 l 

“ b 2 1^7 +2(2b 2 2 p- I -b 1 ) 


+ g ^2iP~ 4 b 2 2e2 ( b 2 2 P~ 1 “ bj ) [ 3 (l~e 2 ) 1/2 - {l~e 2 )~ 1/2 e 2 j 


and 


^24 

<?q, 


24 


Ub; 1 


c9b 2 

Jq7 


- 5 p“ 


.ip 

dq 


;)*A 




[ 4 {l-e 2 ) 1/2 - (l-e 2 )- 1/2 e 2 ] 


23 


By the corrected versions of equations (116.3) and (116.4) of Reference 8, we have: 


B u = - 2PQ + g c i 


13 


24 Ci Q 3 


1 

64 ^2 


C. Q 4 


and 


16 C 1 c 2 0 3 


' 256' C 2 2 


Thus, for i = 1, 2, 3, 


6 B n / 9 , \ <90 (9P 3 

S^7 = (8C,Q 2 -2P)^ - 2Q^ + g ^ 


6q~ 


5B.3 

5q, 


1 / 6C, dQ 

24 Q2 \ Q 5q" +3C i dq~ 


<9B, 


<9C„ 


dQ 


di\ = 64 Q <96, +4C 2 <9q; j ’ 


<9B 


23 


<9C, 


d( h 


( SC 

c 2 q^ +c i Q Sr + 3C i c 2 667 


<?Q 


and 


5B 24 3 , 1 dC 2 dQ 

'5qi " 128 C 2 Q dqj + X 2 <9q. 


TIME-VARYING PARTIAL DERIVATIVES IN THE DIFFERENTIAL CORRECTION 

We shall continue the use of the generalized notation introduced in the preceding section, 
whereby we let q,(i = 1, 2, 3) refer to the orbital elements a, e, and S, respectively. The 
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time-dependent parameters will involve, additionally, the orbital elements /?,, /3 2 , and /? 3 , which we 
shall represent as /3, (i = 1, 2, 3) in the generalized form of the partial derivatives. 

By combining equations (121.1)* and (123.1), we have: 

M s = 2nv x (t +/ 3 , -c 2 / 3 2 a " 1 B/ B/ 1 ) . 

Thus, for i = 1, 2, 3, 

<?M s a 

(2^,) 


- 2 7 TV, C 2 ySj a 2 1 B 2 1 


aB l' , <?a 2 . 5B 2 \ 

dq, B 1 a 2 dq, ' B 1 B 2 dqj ' 


Also: 


and 


^M s 

w x = ■ 

(9M s 

^ - - 277V, c 2 a 2 - 1 B 1 'B 2 -‘ , 


<9M 

^ = 0 • 


In what follows, whenever a partial derivative with respect to is zero, it will not be given. 
By combining equations (121.2) and (123.2), we have: 


\ = 2™ 2 [t+^ 1 +/3 2 a 2 - 1 A 2 - 1 (a+b, +A,)] 


Thus, for i = 1, 2, 3, 

5q7 = [t + /?, +P 2 a' 1 A' 1 (a + b, +A,)] ^ [2nv 2 ) 


+ 2t w 2 P 2 a ~ 1 A 2 1 


5b, 

S H + 


d/± 1 I da 2 c9 A ^ V 

cfqT ~ ( a + b i + A i ) ^“2 * dq 7 + A 2 * dqT" j ' 


*A11 equation numbers, used in specifying the defining relation for a given variable, will henceforth refer to Reference 8, 
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Also: 


W[ = ^ 


and 


dV s 

jf 2nv 2 a' 1 A; 1 (a + b, + A,) 


By equation (131.1), letting £ = M s +E 0 , we have: 


Thus, for i = 1, 2, 3, 


in P Q M. 


dS , O ( ■ p de' 3M S 

^7 " (l-e cosE) + a q . 


Also, for i = 1, 2, 


d£ *M s 

= (l-e' cosP)- 1 ^ 


By the anomaly connections given in equations (132), taken to zeroth order, 


v ' = (cos £-e)(l-e cos £) 1 


and 


;inv' = (l - e 2 ) 1/2 ( 1 - e cos £) 1 


where 


V 0 = v' - M S 


Then, for i = 1, 2, 3, 

<9vq 

d *i 


(l-e 2 ) sin £ + S 2i sin 2 £ ( s in v ' )" 1 ( 1 " e cos 8) 2 ^q. 
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Also, for i = 1, 2, 

dv o I Jg 5M 

W[ = (l-e 2 ) sine(sinv')- 1 (1-ecosSy* jjj- - 

By equation (133), we have: 

V o = “2 u-^A 2 B 2 -‘v 0 . 

Thus, for i = 1, 2, 3, 


dV n 


da 


dq i ( 2a l) 1/2 U 1/2 B 2 1 1^ A 2 v 0 dq i 2 a 2 a l~ 1 A 2 v 0 <9q. " 2 a 2 u 1 A 2 v 0 7^7 


da 


1 1 


<?A 2 _ dB 2 dv Q \ 

a 2 v 0 dqT “ a 2 A 2 B 2 1 v 0 ~dcj~ + a 2 A 2 7>q“ j 


Also, for i = 1, 2, 


"o , . 


W[ = a 2 (-2ai)- 1/2 u' ,/2 a 2 B 2 -' • 

By equation (136.2), we have: 

M l = - ( a+b i) _1 [( A i + c2 A 2 B i' B 2 -1 ) v o + c2a 2* 1 (-2 a i) I/2 u 1/2 B 12 Sin2(f, +»„)] . 

Thus, for i = 1, 2, 3, 

dM, f ( 5b. \ 5v 0 r*A, / 5A. 

5 q 7 = “ ( a+b i) _1 | M i ^ S li + 5 ^~) + ( A i +c 2a 2 b ,' b 2 - 1 ) 75- + v 0 |^7 + c 2 B 2 -‘ (b/ ^7 


5b; 

+ ^2 5q. A 2 ® 


5B 2 V 

*' 82 1 *«>i ) 


+ c J a~ -2a, 


1 . * a i da 2 \ /5V d x V n \] 

2 a i' lB i2 5^ -S 1B 12 (- 2a ,) 1/2 u 1/2 B 12 cos 2(V S + ? 0 )^ + 


Also, for i = 1, 2, 

<?M, 


^ M i ^v 0 /dm dm n \ 

W, = " ( a+b i)' 1 |_( A l + c2 A 2 B l' B 2 _1 ) W, +2c 2 a 2 - 1 (- 2 ai) 1/ 2 u 1/2 Bi 2 cos 2^ + ¥„ ) ^ + jfi-J 
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By equation (137), neglecting terms of third order, we have: 


(l 


£)" 


1 Mj - 


s' (l-< 


;S)- 3 M. 


Thus, for i = 1, 2, 3, 


<5q, 


<?M, 

(1 - e' cos S)‘ 1 “ Mj (1 - e' cos £) 2 

3 

- 2 ~M 1 e / ( 1 - e ' cos 8) 1 


f' 


o ^ 

„ de'\ 


1 

sin8 d^l 

- cosE d^j 

_ Mj ( 1 _ e ' cos £) 3 sin£ 

2 M 

n d£ 

a de'\ 

1 - o dB 

c 9 Mj 

SinE d^~ 

cost dq j 

+ 2 Mj e cot t dq ^ + e 

d % 


Also, for i = 1, 2, 


3E <?M 

- (1 - e 1 cos 8y 1 [l - Mj e ' ( 1 - e ' cos 8) 2 sin 8j 


M x e 


^ p r i 3 i 

( 1 - e ' cos 8)“ 2 Is in 8 + y Mj ( 1 - e' cos 8)' 1 cos 8 - ^ Mj e' (1 - e' cos 8)" 2 sin 2 8j 


By the anomaly connections, taken to first order, 

cos v" = jcos (8 + E j) - ej jl - e cos (8 + EjJ 


and 


sin v" - (l - e 2 ) 1/2 [\ - e cos (8 + EjJ sin (8 + Ej 


where 


Vj - v" - v' - v" - (M s + V 0 ) 


Then, for i = 1, 2, 3, 



+ 


8 2i sin 2 (8 + 


( sin v 



e cos (8 + EjJ 


d v ' 
dq. 


de' 
l (9q x 
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Also, for i = 1, 2, 


dfi\ - (l -e 2 ) sin (S <E,) (sin v")’ 1 [l - e cos (8 + Zj J 2 (j^ + 


In the above and in the following, we employ the symbolic notation: 


dv' dU s dV 0 


<?q 4 dq i 


dv' _ dv o 

W, = BP, + dp. 


By the corrected version of equation (139), we have: 


V l = - B 2 ' [ B 22 Sln 2 (». + *o) + B 21 COS (¥, + »„)] 


(”2a 1 ) 1/2 u 17 2 B 2 1 (A 2 v x + A 2 j sin v ' + A 22 sin 2v') 


Thus, for i = 1, 2, 3, 


r i dB i \ dB o 

B 2 2 [B„ sin 2 (».+»,) +B 21 cos (¥.+? 0 )J ~ B a -W»in2(f.+f 0 ) ^ 


dB 21 fdV dW Q 

+ COS (f, + ¥„) -5^- + [2B 22 cos2(f, +»„) -B 21 sin (< V s + ¥„)] + ^7 


+ (-2aj)‘ 1/2 u” 1/2 B~ 1 


^•2 a 2 +u '‘^ + 2B 2 -ij-| (A 2 Vj 


A 21 sin v' + A 22 sin 2v') + a 2 (“2a J ) ' 1/2 u* 1/2 B 2 1 


' dA 2 
v i 


2 dv i . , (3A 21 

7 + A 2 3^7 +sinv 


■ ^22 . * <9v' 

+ sin 2v , 9 q_ +(A 21 cosv +2A 22 cos2v'J 
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Also, for i = 1, 2, 


/5T s 

b; 1 [2B 22 cos2(¥,+¥ 0 ) -B 21 sin (*.+¥ 0 )J ^ ) 


r ^ v i \ <9v'i 

( _2a i)" 1/2 u " 1/2 B 2~ 1 A 2 W~ + ( A 2 1 cos v ' + 2A 22 cos 2v') 


The following time -dependent partial derivatives are used only when the differential correc- 
tion includes periodic terms through the second order. 

By equation (143.2), we have: 


( a + b i)~ 1 { A i v i + A n si 


sin v ' + A 12 sin 2v' 


: 2 a 2 1 (-2a I )‘ /2 u 1/2 [b/ 4*, +B,, cos (q» s + ¥ 0 ) 

+ 2B 12 'l- 1 cos2('P s+ 'l« 0 ) + B 13 cos 3('1< s + >P 0 ) + B 14 sin 4 (v s + ¥ 0 )]} 


Thus, for i = 1, 2, 3, 


c?M 2 / db\ (9Aj dv t 1 1 ^A 12 


( a + b 1 <* 2 S h + + A i ^ + sinv '-^7 + sin 2 v # 


(A 11 cosv'+2A 12 cos2v') + c 2 a 2 " 1 (~2a t ) 1/2 u 1/2 ^ u" 1 3^7 

1 da i da 2 \ r 

2 - -2 1 L B i ,,, i + »„«»(¥.+¥„) + 2 B 12 cos 2 (¥ s + ¥ 0 ) 


B 13 cos 3 (S’ s + +B 14 sin4(«P s + ¥ 0 )] + c 2 a; 2 (-2a,) 2/2 u " 2 U 


dV <9B I dB dV\ 

+ B i 5^7 + cos (¥,+¥„) + 2 K^“ + B 12 ^-J cos 2 ('P s + ¥„) 


<7 0 J 3 (7 Jt> 14 - 

+ cos3(¥.+¥ 0 ) -55- + sin 4 (¥ s + ¥„) - [b, , sin (¥. + ¥ 0 ) 


4B 12 'P 1 sin2(4- s+ <P 0 ) + 3 B 13 sin 3(f. + ¥ 0 ) - 4 B 14 cos 4(>P s + ¥ 0 )J 


a<P s <9>P 0 
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Also, for i = 1, 2, 


5M 2 

dK 


dv 


( a + b j ) ” 1 U, a/. + (A n co S v' +2A 12 cos2v / ) 


c 2 a 2 1 (-2aJ 1/2 u 


1/2 


B i 5/3. +2B 12 cos2(f, +»„) - [Bjjsin (*, + ¥„) 


4B i^> sin 2 (*P s + >P 0 ) + ^B 13 sin 3 (>P s + *P 0 ) - 4 B„ cos 4(¥, + ¥„)] ^ 


By equation (143.1), we have: 


E 2 _ j^l - e 1 cos (s + Ej)J 1 1 


Thus, for i = 1, 2, 3, 


dE 2 _ r . 1 3M, 

dq . = [l-e' cos (e + EjJ ^7 


- M. 


^1 - e ' cos (8 + Ej)j 


( 8+E l) (/q^ + 5^) "cos (S + Ej) f§- 


Also, for i = 1, 2, 


^ |T ^ j^j 

jf, [l-e' cos (S +El )] _1 ^ -M 2 [l-e'cos (8 + Ej )] " 2 e ' sin (g + £,) • 


Since E = 8 + Ej +E 2 , then, for i = 1, 2, 3, 


Also, for i = 1, 2, 


5E _ _dS_ 

" dq i + <9q i + 5q“ 


5E _ _dS_ £^i ^£2 

5/3. = 5/S, + 5/3. + 5/3. ' 


By the anomaly connections, taken to second order (where we choose the notation v rather 
than v "'), 


iv - (cosE-e) ( 1 -ecosE)' 


31 



and 


sin v = (l - e 2 ) 12 (l-ecosE) *sinE 


where 


V - v" - V - (M. +v 0 + v,) 


Then, for i = 1, 2, 3, 


<5q ( 


<?E 


(l - e 2 ) sinE + S 2; sin 2 E (sinv)”*(l-e cos E) 2 


+ d %) 


Also, for i = 1, 2, 


dv „ 


J7/~ - (l - e 2 ) sinE(sin v)' 1 (1 - e cos E) 2 


d E fdv' dv i\ 

W[ ~ ^ + afij 


Since v = v' +v t +v 2 , then, for i = 1, 2, 3, 


dv _ dv ' dv i d y 2 
’5q‘ “ dq i + dq“ + dqj 


Also, for i = 1, 2, 


dv <9v' dV l av 2 

d/3. - d/3. + d/3. + d/3. 


By the corrected version of equation (145), we have: 

*a = - B 2 _1 [- B 2i*i sin (V. + ’o) +2 B 22 'P 1 cos2('P s+ 'P 0 ) +B2 3 cos3(?,+?o) 

+ B 24 sin 4 (v s + ¥ 0 )J + a 2 1/2 u" 1/2 B 2 “ 1 ( A 2 v 2 + A 2J v x cos v' 


+ 2 A 22 v i cos 2v' + A 23 sin 3v' + A 24 


sin 4v ') 
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Thus, for i = 1, 2, 3, 


3qf = B 2 _ 2 [ _B 21 , 'l Sin (^s +,, o) + 2B 22 , 'l COS 2 (’ P s + \) +B 23 COS 3 (' i 's +,P o) 


<9<P, 
21 <?q i 


(*. + *o) 


T <JB, f / 5B-, 

B 24 s * n 4(^5 + *^ 0 )J 

( dB 22 d'VX dB $B 

2 ^i ~d^~ +B 22 d^~) cos 2(<P s + -P 0 ) + cos3(? s +w 0 ) - 5 ^- +Sin4(<p s +w 0 ) 


[ B ai*i cos C*. + *o) +4B 22 ,p i sin2(>P s +¥„) +3B 23 sin 3 (?. + >P 0 ) 




4B 24 cos4(<Ps^ 0 )J (”2 a j ) 1/2 u" 1/2 B 2 ~ 1 


da 


2 1 f , 3<l 1 , <?U 

5q. - 2 a 2 ( a l dq. + u dq. 


dB ? 

+ 2 B 2 * dq. 


( A 2 v 2 +A 2 i v l cos v ' + 2A 22 v i cos 2v ' + A 23 sin 3v' + A 24 sin 4v ') 


a 2 (*2a 1 j“ 1/2 u 1/2 B 2 1 


<?A. <?v. 


aA„ 


2 w v 2 / un 21 ^ V l\ 

v « + A 0 XT" + (v, “XT" + A,, ~Sq~J 


2 dq. T A 2 d q. + l v i dq. + A 2 1 ^T/ cos v 


+ 2 v 


«9 A 22 _ <3 Vj 


dA„ 


dA 


l dq, + A 22 di~) cos 2v ' + sin 3v' + sin4v' ~g^ 


( A 2 1 v i S ^ n v< +4 A 22 v i s * n 2v' 3 A 23 cos 3v' - 4A J4 cos 4v') r 


Also, for i = 1, 2, 


dv 2 I r av 

= B 2'M[ B 21 sin K+’O - 2B 2 2 COs2 ('P s + 'Po)] W: + [ B 2l’ P l COS (^ , + » 0 ) 


(9^ \ 

4B 22 'P 1 sin2(«P s + «P 0 ) +3B 23 sin 3 (*. +? 0 ) - 4 B 24 cos 4(«P s + *„)] ( ^ , 


r <?v d v 

a 2 (-2^i) _1/2 1/2 B - 1 | A 2 + (A 2 , cos v* + 2 A 22 cos 2v') jf; 


- (A 21 v, sin v' + 4 A 22 v, sin 2v' - 3 A 23 cos 3v' - 4 A J4 cos 4v') 


dv ' 
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Since W + + ^1 + ^2 > ^ en ^ or 1 ~ 2, 3, 


<94* 


(9W s aW 0 aw x <?W 2 

(9q x + dq i + + ^ 


Also, for i = 1, 2, 


^ ^0 ^1 ^2 

a^ = a/3. + a/3. + a/3. + a/3. 


This concludes the evaluation of the partial derivatives of the uniformizing variables E, v, 
and W when the computation is followed through terms of the second order. If, however, second- 
order precision is not required, then the partial derivatives of M 2 , E 2 , v 2 , and w 2 may be omitted, 
and the above partial derivatives of the uniformizing variables reduce to 


and 


<?E 

dF 

**1 


del, 

+ 


<?v 

<?v ' 

dv i 


aq. 

+ 


aw 

d % 


aw aw n aw. 

3qf + 


Likewise for the partial derivatives of E, v, and w with respect to ( 3 . (i = 1, 2). 

We now continue with the necessary equations preparatory to determining the partial deriva- 
tives of the inertial co-ordinates X, Y, and Z. 

By special cases of equations (104) with y = W± (l/2)n 9 we find: 


and 


cos E 2 ' ^ (l-e 2 sinW) 1 (e 2 - sin'P) 

sinE 2 ' = (l~e 2 2 j l/2 (l - e 2 s in w) “ 1 cos W , 

cos Ej = (l+e 3 sinWj _1 |e 3 + sin , 
sinE 3 - _ (l _e 3 2 ) 1/2 (l + e sin w) ~ 1 cos W 
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Thus, for i = 1, 2, 3, 




dq. " (1- e^in'P)' 1 


I, 2 x 1/2 d¥ . ,,-1/2 

( X ~ e 2 ) - (* _ e 2 ) COS 


ae 2 ] 


and 


dE± 

d % 


(l + e 3 sin¥)-‘ (1 -e 3 2 ) 1/2 t (l - e/) 


- 1/2 d e 3 

cos,p 


Also, for i = 1, 2, 


**2 


d/3. = (1 - e 2 sirtV)- 1 (1 - e 2 2 ) 1/2 


and 


dE 3 ' 


d/3. = (l +e 3 sin'P)' 1 (1 -e 3 J )‘ 2 


By equations (114.1) and (114.2), we have: 

*0 = \ (1 - 2C)' l/2 E 2 ' + \ (1 + 20" 1/2 E 3 * 

and 

*i = | (l - 20' 1/2 e 2 ' - § (1 + 20" 1/2 e 3 ' • 


Thus, for i = 1, 2, 3, 


dy 0 a dE-' i dE/ i r ^ r 

d^ = 2 (1-2;)‘ I/2 ^- + 2 d+20" 1/2 + 2 [E 2 *(l-2Sr 3 ' 2 - E 3 '(1 + 2C)- 1/2 ]^ 


and 


d Xi i d E J i d E J i r ^ 

d^“ = 2 (1-20‘ 1/2 ^ - | (1 + 2S)' 1/2 -^- [e 2 ' (1 -2C)' 3/2 + E 3 '(l+20' 3/2 ]-^7 
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Also, for i = 1, 2, 


d *o 1 


af. = j (1 - 20" 1/2 ~sf. + f (1 + ’ 


dE ' 


2 1 


dE,' 


and 


1 ^E' 1 <9E 3 

2 (f " 2 0" 1/2 ~d]3. - 2 d + 2 0’ 1/2 ^- 


<?Xi 

By equation (150), we have: 


4 > ~ /3 3 ~c 2 a 3 ( _ 2a : ) 1/2 (A 3 v + A 31 sin v + A 32 sin 2v+A 33 sin3v + A 34 sin 4v) 

+ a 3 a 2 -!u 1/2 |(1 -S)“ l/2 [(h, + h 2 ) * 0 + ( h !- h 2 ) Xj] 

+ B 3 V- | CjCjQcos'P + ^ C 2 2 Q 2 sin 24>j 

Thus, for i = 1, 2, 3, 


_ 2 
d q : " ~ c 


da. 


da. 


( _ 2 a j) 1/2 5q~ + a 3 (”2 a i) * 3/2 

A 3 3 sin 3v + A 34 sin 4v) - c 2 a 3 ^-2a J j 


(A 3 v + A 3 j sin v + A 32 sin 2v 


1/2 


dA 3 dA 31 dA 3 2 

v dV + sinv d^T + sin 2v 


dq. 


dA , 


dA 


33 34 

+ s in 3v + sin 4v ^ + (A 3 + A 3 x cos v + 2 A 3 2 cos 2v + 3 A 33 cos 3v 


v dv 

4A 34 cos4v) 


+ a, 1 u 


l n i/2 


da, 

Jq 


f ' a 3 a 2 1 + T a 3 U ~ 1 {( 1 " S >" 1/2 [(h, + h 2 ) 


X 0 


( h i _h 2 ) Xj] + B 3 'J<- f Cj C 2 Q cos V + ^ C 2 2 Q 2 Sin2>p| 
a 3 a 2 - 1 u^ 2 j~ 8 3i (1-S)-^ 2 [(h 1+ h 2 ) X 0 + ( h l' h 2 ) X,] + O ~ s ) 


( h i +h 2 ) 


5 X 0 

<?q ; 


d *l 


ah, 


dh„ 


( h i _ h a) dq] + (x 0 + Xi) dq] + (x 0 " X,) d q * 
dC„ 


dB, 


3 3 / “-i 

+ ^ 5 q 4 ( ^ 2 Q dq. 


dC, 


<?Q\ 3 „ „ dc 2 . „ dQ 


+ C i Q ^7 +C 1 C 2 d^J co ^ + l6 C 2 Q(Q^7 + C 2 a^i sin 


2<P 


/ 3 3 \ <9>p] 

( B 3 + 4 Cl C 2 Qsin-P +T g C 2 2 Q 2 cos 2 pJ 
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Also, for i = 1, 2, 

W[ = " c2a 3 (' 2a i) 1/2 ( A 3 +A 31 COSV+2A 32 COs2v+3A 33 COS 3v + 4A 34 COs4v) jj 


+ a. a, 1 u 1/2 2(1 


f, dX 0 dXi 

-S)-i ^(h, + h 2 ) + (h^h,) 


( 3 3 \ 

B 3 + TC 1 C 2 Qsin«l< + jgC 2 2 Q 2 cos2 v) Jp~, 


and 


d<t> 

a/s 3 = 1 • 


Since the spheroidal co-ordinates are given by 


p - a (l - e cos E) 


and 


r] - P + Q sin , 


we then have 


Bjp 

da 


1 


e (cos E - a sin E 


dp _ / dE \ 

- a ^esinE^ - cosE J , 


and 


ap . . „aE 

as aesinE as ‘ 


Also, for i = 1, 2, 


dp . „ dE 

W = aesinE a/T • 
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Further, for i = 1, 2, 3, 


dp dP dQ dV 

d^ - + * in *5^7 + 


and, for i = 1, 2, 


-TL = Q cos»P — 


d l _ 

dP i 




The partial derivatives of p and v with respect to /3 3 are both zero. 
Finally, the inertial rectangular co-ordinates are given by 


X - }/{p 2 + C 2 ) (l - Tj 2 ) COS <£ 


Y - }/{p 2 + c 2 ) (l - 7] 2 ) sin ^ 


and 


Z = P7) - 5 


Thus, for i - 1, 2, 3, 


dX 


dq. 


= X 


ay 


dq. 


['I''*' 1 ) '557] 


and 


dZ _ drj dp 

dq” " p dq~ + 77 dq“ 


Also, for i = 1, 2, 



and 


and 



dZ _ dr} dp 

d/3. ~ P d/3. + V d(3. 


dX 


-Y , 


d Y 

w, - x • 


az 

~nr - o 


OBSERVATIONS OF THE RIGHT ASCENSION -DECLINATION TYPE 

The differential coefficients in the form of partial derivatives of the inertial rectangular co- 
ordinates with respect to the orbital elements are completely general in the sense that they are 
functions only of the mathematical theory of orbital satellite motion. Thus, they are applicable to 
observations of spacecraft position recorded in any format whatsoever. Previously, we had as- 
sumed that the observational data were recorded as direction cosines with respect to a topocentric 
latitude-longitude-zenith co-ordinate system. Another format frequently used for recording ob- 
servational data is the right ascension-declination type. In this section, we shall discuss the minor 
variations in the equations that arise when this type of data is utilized. 

The co-ordinate system adopted for the use of right ascension -declination data is also situated 
at the tracking station on the Earth’s surface, but its three co-ordinate axes are parallel to the 
respective axes of the inertial system. Again designating the topocentric ’’local 11 co-ordinates by 
the subscript "M", we have, in this case, the Z M -axis parallel to the Earth’s polar axis, and the 
\ -y m plane parallel to the equatorial plane of the Earth. The x M ~axis extends toward the vernal 
equinox, with the Y M -axis orthogonally to the east to form a right-handed system. 

The observed right ascension, a 09 is measured eastward from the vernal equinox, and the ob- 
served declination, 8 0 , is measured as positive above and as negative below the equatorial plane. 
The corresponding computed values of the right ascension and the declination are given in terms 
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of the local co-ordinates by 


a 

c 


arctan 



and 


8 


arctan 


(V + V) 


1/2 


In order to obtain a satellited local co-ordinates from its inertial co-ordinates, the inertial 
co-ordinates of the observation point at the time of observation, denoted (x T , y t , z t ), must be 
known. However, no rotations are necessary to bring the two systems into coincidence in this 
case, since the topocentric and inertial co-ordinate systems are parallel. A single translation 
will suffice. Hence, the relations for the local co-ordinates of the satellite are simply 

X M = x - x T , 
y M = y - y t , 


and 


z m Q z ‘ Z T • 

Notice that the above simplified relations are obtained from those of the direction-cosine-data 
case by the artifice of setting «A X = 0 and 0 D = tt/ 2 in the corresponding equations for x M , y m , and 
Z M given earlier. 

The first-order Taylor’s series expansion for the equations of condition corresponding to each 
time of observation are 


and 



i = 1 
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where q A ( i = 1, 2, • • • , 6) are the orbital elements. Expanding the above partial derivatives by 
the chain rule yields 


and 


da da 3X U da d Y u d a d Z u 

^7 r 


a qi dX^ dq. + <9Y M dq, + d Z„ d q ; 


From the expressions for a c and S c in terms of the local co-ordinates, we find 


da 

dT u = 

r M 2 

d a 

= 


d a 

0 , 

33 c 

* x m 

XmZm 

" 

38 c 

Y m Z m 

* Y « 

r M R M 2 


and 

d 8 c _[m 

' + V ' 

where R M - (x M 2 + y m 2 + zj ) x/1 and r M - (V + V)^- 

Since the station co-ordinates, X T , Y T , and Z T , are independent of the orbital elements (and 
merely geodesic functions), the following simple relations apply: 


d Qi 


ax 
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and 


dq i 


dY 

d *i 


dZ M JZ_ 
dq i ' (9q i 


The differential coefficients, dX/dq i9 dY/dq i , and dZ/dq i (i = 1, 2, • • • , 6) are precisely those 
that have been evaluated previously in the differential correction scheme. 


REMARKS 

The differential correction process removes inaccuracies of the initial conditions (the nominal 
observations) and accounts for the effects of forces not considered by the analytical orbital theory. 
Such neglected forces may include aerodynamic drag, electromagnetic effects, solar radiation 
pressure, meteoric bombardment, and residual gravitational influences (including those arising 
from lack of spherical symmetry in the satellite, as well as perturbing planetary potentials). This 
is all accomplished by producing a mean set of orbital elements through an iterated least-squares 
fitting of the first-order Taylor T s series expansion of the conditional equations to numerous ob- 
servational values. Generally speaking, the fitting will improve as greater numbers of observations 
are considered and as the time span represented by the observational data is lengthened. However, 
the complexity of the mathematical processes involved in the fitting increases rapidly as additional 
observations are admitted. Because of this latter constraint, it is often advisable to perform re- 
peated differential corrections at various intervals of time (known as "epochs") rather than at- 
tempt to accommodate all the data in a single iterated fitting. 

The orbital improvement method of differential correction discussed in this paper has been 
programmed, primarily in the FORTRAN language, for use on a large-scale electronic digital 
computer. The analytical nature of the entire procedure assures a very rapid computational 
process. Application of the differential correction (combined with an orbit generator of position 
and velocity components) has been made to both actual observational data of artificial Earth sat- 
ellites and to artificially generated "data" for extremum cases, e.g., polar and equatorial orbits. 
The results have proven entirely favorable. Particular experimental applications will be published 
in a later paper. 
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APPENDIX 


FUNCTIONAL DEPENDENCIES OF THE 
PARTIAL DERIVATIVES 


The following tables are intended to display the inter-relationship of the various parameters 
whose partial derivatives appear in the differential correction. The partial derivative of each 
parameter in the left column is seen to be a function of those, and only those, partial derivatives of 
parameters occurring in the respective line of the right column. There is a certain amount of flex- 
ibility in the ordering of the partial derivatives occurring in the differential correction, as is 
demonstrated by the functional dependencies illustrated by these tables. 

Table I 


Time-independent Partial Derivatives 
(All taken with respect to orbital elements a, e, and S.) 


Partial Derivative 

Functional Dependence on Other 
Partial Derivatives 

p 

(none) 

b i 

p 

b 2 

P, b j 

P 0 

P, bj, b 2 

a 2 

Po 

( a o Po)' 1 

b i- Po 

u 

Po ' ( a o Po)' 1 

C 2 

( a o Po)' 1 - u 

C ! 

Po - U - C 2 

P 

P0’ U ’ C 2 

“l 

b 1 

/ b l\ 

\ b J 

b l- b 2 

(V\ 

\ p / 

P, b 2 

A , 

Ab\ /bA 

P ’ \ b 2 / \P/ 
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Table I (Concluded) 


Partial Derivative 

Functional Dependence on Other 
Partial Derivatives 

*A 

24 

P. b 2 

* B n 

C,, P, Q 

* B ,3 

C i> Q 

* B 1 4 

c 2 ,q 

* B 2 3 

C 2 , Cj, Q 

* B 2 4 

c 2 , Q 


Table n 

Time-varying Partial Derivatives (All taken with respect to orbital elements 
a > e > s > and f J) 2 - Exception: i- , X, and Y are taken with respect to d 3 , as well.) 


Partial Derivative 

Functional Dependence on Other 
Partial Derivatives 

M 

s 

V B/, B 2 , (2^) 


b r 2 ' A r A r C 27 '^) 

F 

e\ M 

s 

v o 

M , F 

s 

'J'o 

a 2’ U ’ V A 2* B 2 ’ v 0 

M 1 

b,, a 2 , u, a,, Aj, A 2 , Bj'. B 2> B 12 , «|< s , v Q , 

E ! 

e ' , c, M , 

V ' 

M s , v 0 

V 1 

8, Ej, v' 


a 2 , u, a Jt A 2 , A 21 , A 22 , B 2 . B 2 ,, B 22 , » 0 , v', v, 

*m 2 

b r i J1 u,a 1 ,A I1 B',B lll A 11 ,A 1J .B 111 B 131 B I4 . 


f,. 4<o- v '- v i- 

*e 2 

e', e, Ej, M 2 

E 

e, e,,e 2 

*v 2 

v' , V r E 

V 

v', V,, V 2 
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Table II (Concluded) 



Note: In Tables I and 13, the asterisk indicates partial derivatives of 
parameters used only if the differential correction includes 
periodic terms through the second order. 




